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We have studied the effects of the 2006 El Nino on tropospheric O 3 and CO at 
tropical and sub-tropical latitudes measured from the OMI and MLS instruments on the Aura 
satellite. The 2006 El Nino-induced drought allowed forest fires set to clear land to bum out of 
control during October and November in the Indonesian region. The effects of these fires are 
clearly seen in the enhancement of CO concentration measured from the MLS instrument. We 
have used a global model of atmospheric chemistry and transport (GMI CTM) to quantify the 
relative importance of biomass burning and large scale transport in producing observed changes 
in tropospheric O 3 and CO. The model results show that during October and November both 
biomass burning and meteorological changes contributed almost equally to the observed increase 
in tropospheric O 3 in the Indonesian region. The biomass component was 4-6 DU but it was 
limited to the Indonesian region where the fires were most intense. The dynamical component 
was 4-8 DU but it covered a much larger area in the Indian Ocean extending from South East 
Asia in the north to western Australia in the south. By December 2006, the effect of biomass 
burning was reduced to zero and the observed changes in tropospheric O 3 were mostly due to 
dynamical effects. The model results show an increase of 2-3% in the global burden of 
tropospheric ozone. In comparison, the global burden of CO increased by 8-12%. 
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1 Introduction 


El Nino and La Nina events are major sources of inter-annual and decadal variability in tropical 
tropospheric ozone (O 3 ) [e.g., Ziemke et al. 2003]. The atmospheric effects of El Nino events are 
generally a change in convection in the tropical troposphere associated with an eastward shift of 
the warm SST anomaly and large scale Walker circulation. This shift results in an increase in 
tropospheric column O 3 in the western Pacific and a decrease in the eastern Pacific relative to 
non-El Nino years. The effect of El Nino on specific humidity is usually opposite to that of 
column ozone [e.g., Chandra et al., 1998; 2007], During La Nina years, dynamical processes are 
largely reversed w hich results in a decrease (increase) of column ozone (specific humidity) in the 
western Pacific. The effects of El Nino on tropospheric composition have been extensively 
studied from both satellite and ground based measurements [e.g., Chandra et al., 1998; Fujiwara 
et al., 1999; Thompson et al., 2001] and by using global models of atmospheric chemistry and 
transport [e.g., Peters et al., 2001; Sudo and Takahashi, 2001; Chandra et al., 2002, 2003; 
Duncan et al., 2003; Zeng and Pyle, 2005; Doherty et al., 2006] 

Recently Logan et al. [2008] studied the effects of the 2004 and 2006 El Nino events on 
tropospheric profiles of CO, O 3 , and H 2 O measured from the Tropospheric Emission 
Spectrometer (TES) on the Aura satellite. Their findings were generally consistent with the 
observed characteristics of O 3 and FLO inferred from previous El Nino events. The changes 
during the 2004 El Nino in tropospheric O 3 and LLO from TES are similar to those inferred from 
Ozone Monitoring Instrument (OMI) and Microwave Limb Spectrometer (MLS) flown on the 
same satellite [ Chandra et al., 2007], The 2004 El Nino was a weak event compared to 2006 
and significantly weaker compared to the El Nino of 1997-1998. This was reflected in the TCO 
(tropospheric column ozone) anomaly over the Indonesian region in the western Pacific as 
shown in Table 1 of Logan et al. [2008]. For example, the TCO anomaly in November 1997 
inferred from EP TOMS was 14.4 Dobson Units (DU; 1 DU = 2.69 x 10 20 molecules-m 2 ). The 
corresponding values from TES for November 2004 and 2006 were 6.6 DU and 11.1 DU. Logan 
et al. [2008] attributed most of these changes to differences in the magnitude of CO emissions 
from fires in Indonesia which in these years were estimated to be 193 Tg, 24 Tg, and 82 Tg, 
respectively. They also suggested significant contribution from NO x production due to lightning 
in late November and December 2006 when CO production due to large scale forest fires 
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decreased significantly. The effects of forest fires in the Indonesian region during the 1997 El 
Nino was studied by Chandra et al. [2002] using the GEOS-CHEM global model of chemistry 
and transport. Their study suggested that about half of the increase in tropospheric column 
ozone in the Indonesian region was due to biomass burning and the other half due to dynamical 
effects. A similar conclusion was arrived by Sudo and Takahas hi [2001] from their model study. 
Their conclusion about the contribution of biomass burning was indirect since they did not 
explicitly include the effect of biomass burning in their model study. 

The combined OMI and MLS instruments on the Aura satellite have been providing near global 
measurements of TCO from August 2004 to present [e.g., Ziemke et al. , 2006; Schoeherl et al., 

2007] as discussed in Section 2. In addition the MLS instrument has been providing daily 
measurements of CO at several pressure levels in the troposphere and stratosphere [Livesey et al., 

2008] . The purpose of this paper is to study the effects of the 2006 El Nino on tropospheric 
composition measured from the OMI and MLS instruments on the Aura satellite and evaluate 
specific roles of the various processes using a global model of chemistry and transport. The 
global model used in this study is the NASA GMI CTM (Global Modeling Initiative Chemical 
Transport Model). It successfully simulates a wide range of observations of chemical 
constituents in the troposphere and stratosphere, including data collected by instruments on 
board satellites. We provide a brief overview of the relevant tropospheric processes in the GMI 
CTM in section 3. In our study, the El Nino related changes in O3 and CO fields based on OMI 
and MLS measurements are compared with corresponding changes based on GMI CTM 
simulation to delineate the relative importanee of biomass burning and large scale transport. 

In the following, section 2 discusses the satellite data, section 3 the GMI model, and sections 4 
and 5 compare observed and measured O3 and CO, respectively. Section 6 discusses the impact 
of biomass burning emissions over Indonesia during the 2006 El Nino, and section 7 compares 
the strength and frequency of the 2006 and recent El Nino events with previous El Nino events. 
Finally, section 8 provides a summary. 


2 Measurements from OMI and MLS 
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OMI is one out of a total of four instruments onboard the Aura spacecraft which is flown in a 
sun-synchronous polar orbit at 705 km altitude with a 98.2° inclination. The spacecraft has an 
equatorial crossing time of 1:45 pm (ascending node) with around 98.8 minutes per orbit (14.6 
orbits per day on average). OMI is a nadir-scanning instrument that for visible (350-500 nm) 
and UV wavelength channels (UV-1: 270-314 nm; UV-2: 306-380 nm) detects baekscattered 
solar radiance to measure column ozone with near global coverage (aside from polar night 
latitudes) over the Earth with a resolution of 13 km x 24 km at nadir. Besides ozone OMI also 
measures cloud-top pressure, aerosols and aerosol parameters, NOa, SOa, and other trace 
constituents in the troposphere and stratosphere [Level t et al., 2006]. Total ozone from OMI is 
derived from the TOMS version 8 algorithm. A description of this algorithm may be obtained 
from the TOMS V 8 CD DVD ROM, or from the OMI Algorithm Theoretical Basis Document 
(ATBD) (from the webpage http://toms.gsfc.nasa.gov/version 8 /v 8 toms_atbd.pdf). 

There are two standard OMI total ozone products. Here we use data from the collection 3 
OMT03 version 8.5 (v8.5) product that is based on the Total Ozone Mapping Spectrometer 
(TOMS) version 8 (v 8 ) total O 3 algorithm. A description of this algorithm may be obtained from 
the TOMS v 8 DVD or the OMI Algorithm Theoretical Basis Document (ATBD) from the 
TOMS web page htt p://toms.gsfc.na.sa.gov/version 8 /v 8 toms atbd.pdf . A main difference 
between v 8 and v8.5 is that v8.5 uses measured cloud pressure whereas v 8 and earlier versions 
used an infrared-measured climatology for cloud pressure. Another OMI algorithm based on the 
Differential Optical Absorption Spectroscopy Method (DO AS) technique gives similar estimates 
of total column ozone [Kroon et al., 2008b]. These total column O 3 retrievals have been 
compared with ground-based data [Me Peters et al., 2008] and aircraft-based measurements 
[Kroon et al., 2008a]. 

O3 profiles from the Aura Microwave Limb Sounder (MLS) are used to estimate stratospheric 
column O 3 (SCO). The SCO is subtracted from OMI total column O3 from OMI to yield 
tropospheric column O 3 and mean volume mixing ratio [e.g., Ziemke et al, 2006]. The mean 
volume mixing ratio (units ppbv) is determined by taking tropospheric column O3 (in Dobson 
Units, DU) and dividing this by pressure difference (in units hPa) between surface and 
tropopause, and then multiplying this by the factor 1.27. Validation of MLS \ 2.2 measurements 
is discussed by Froklevam et al. [2008] and Livesey et al. [2008]. The measurements of mean 
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tropospheric O 3 from OMI/MLS and upper tropospheric CO from MLS are studied and 
compared with simulations from a global chemical transport model (discussed next). 

3 NASA Global Modeling Initiative Chemical Transport Model 

The GMI CTM is described and evaluated in Ziemke et al. [2006], Schoeberl et al. [2006], 
Stratum et al. [2007] and Duncan et al. [2007, 2008], The chemical mechanism includes 117 
species, 322 chemical reactions, and 81 photolysis reactions to simulate tropospheric and 
stratospheric chemistry. It includes a detailed description of CVNOx-hydrocarbon chemistry, 
updated with recent experimental data. The chemical mass balance equations are integrated 
using the SMVGEAR II algorithm [Jacobson, 1995]. Photolysis frequencies are computed using 
the Fast-JX radiative transfer algorithm (M. Prather, personal communication). The algorithm 
treats both Rayleigh scattering as well as Mie scattering by clouds and aerosol. The CTM 
simulates the radiative and heterogeneous chemical effects of sulfate, dust, sea-salt, organic 
carbon and black carbon aerosol on tropospheric photochemistry. 

Biogenic emissions of isoprene and monoterpenes are calculated on-line as described in 
Guenther et al. [2006]. The lightning NOx parameterization fits the relationship between model- 
calculated convective mass fluxes and observed cloud-to-ground flash rates (D. Allen and K. 
Pickering, personal communication). Time-appropriate anthropogenic and biomass burning 
emissions include, surface emissions from industry/fossil fuel, biomass burning, biogenic and 
biofuel combustions and contributions from aircraft emissions. They are based on the Global Fire 
Emission Database, version 2 (GFEDv2) and discussed in the GMI CTM papers listed above. 
Monthly total biomass burning of CO and NOx are shown in Table 1 for the later part of 2006 
when the Indonesian fires were most intense. 

Table 1 indicates that the Indonesian fires increased the CO burden by 69 Tg during the three 
month period from September to November 2006. It is about 27% of the global CO burden (251 
Tg) during the same period. The increase in NOx due to Indonesian fires during the same period 
was 1.19 Tg compared to the global burden of about 19.4 Tg. 

The meteorological fields that drive transport are from the Goddard Modeling and Assimilation 
Office (GMAO) GEOS-4 data assimilation system (GEOS-4-DAS) [Bloom et al., 2005]. 
Features of the circulation, such as anticyclones, in the GMI CTM are realistically represented in 
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the simulations because of the data constraints on meteorological analyses. The GEOS-4-DAS 
fields have been regridded to 42 vertical levels with a lid at 0.01 hPa. The horizontal resolution 
is 2.5° x 2° (longitude by latitude). 

A validation of the global model used in this study is given in the Appendix. This analysis 
compares the zonal and seasonal variability of TC'O derived from the model with observed 
variability from OMI/MLS. 

4 Model comparison of El Nino related changes in O 3 with satellite measurements 

As in Chandra et al. [2007] and Logan ct al. [2008], we have chosen October, November and 
December 2005 as a baseline for estimating El Nino-related changes in these months in the 
preceding and the following year. Year 2005 was a neutral year from an El Nino perspective. 
Maps of tropospheric O3 mean mixing ratio are shown in Figure la for October 2004, 2005, and 
2006. They show relative changes in ozone in the El Nino years (upper and lower panels) with 
respect to a non El Nino year (middle panel). All three panels show regions of minimum values 
near the dateline. However, both the upper and the lower panels show a greater eastward shift 
across the dateline which results in an increase of about 6-15 ppbv (4-10 DU) over the western 
Pacific region and a decrease of smaller magnitude east of the dateline. The increase in O3 mean 
mixing ratio is larger in October 2006 (lower panel) which is consistent with a relatively stronger 
El Nino event in 2006 compared to 2004. The model (Figure lb) captures most of the zonal 
characteristics of TCO shown in Figure la, but it tends to underestimate the observed values by 
5-10 ppbv over most of the region between ±10 deg. (Figure lc). Outside this region the sign is 
reversed and the model is biased higher by about the same magnitude. These pattern shifts in the 
model and observation are comparable in November and December months (not shown) and are 
similar to those shown in Figure 2 of Ziemke et al [2000] for the 1997 El Nino event. 

The El Nino related changes in ozone mean mixing ratio in 2006 are clearly discernible in 
Figures 2a, 2b and 2c. These figures show differences in OMI/MLS O3 (top panels) and the 
model O3 (bottom panels) between years 2006 and 2005 for the months of October-December. 
The model and observed values are generally in good agreement except in the western Pacific 
region in October (Figure 2a) where the model values are biased higher. This bias is also 
reflected in CO as indicated in section 5 and is probably due to overestimation of the effect of 
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biomass burning in the Indonesian region in October 2006. In general, observations and model 
indicate excellent agreement in characterizing both broad and small scale features including the 
dipole nature of ozone anomaly in the tropical Pacific region (Figures 2a and 2c). The large 
negative anomaly over eastern Africa/Indian Ocean in December 2006 seen in the data (Figure 
2c) is a manifestation of dynamically induced changes and is well produced in the model. We 
have analyzed the model results to estimate the contributions of the upper (above 500 hPa) and 
the lower troposphere in producing the El Nino related changes in tropospheric ozone. In 
general, the contributions from the upper and the lower troposphere are almost equal. However, 
the large negative anomaly in December in the eastern Africa/Indian Ocean (Figure 2c) comes 
mostly from the upper troposphere. 

5 Model Comparison of CO with satellite measurements 

The El Nino events of 2004 and 2006 caused droughts in the Indonesian region which allowed 
set fires to spread uncontrolled in both October and November of these years. The MLS 
measurements show a relatively larger increase in CO in 2006 as compared to 2004 (figure not 
shown). This is consistent with CO measurements from TES [Logan et al., 2008], The model 
results show CO to be well mixed in the free troposphere even though convection was greatly 
reduced in this region. During October 2006, the model shows an increase of about 200 ppbv in 
the Indonesian region over much of the free troposphere. Figures 3a-3c show differences in 
upper tropospheric CO from MLS (top panels) and the model (bottom panels) between 2006 and 

2005 for October-December. In October, both the model and measurements show elevated 
values of CO in 2006 over the Indonesian region in the western Pacific. One month later in 
November (Figure 3b) upper tropospheric CO from both MLS and the model are greatly 
reduced. They become negligible by December as the fire subsides in this region. CO for both 
measurement and model exhibit very similar spatial and temporal patterns except in October 

2006 (Figure 3 a). The model shows a sharp horizontal gradient with CO mixing ratio exceeding 
200 ppbv in the Indonesian region. The MLS values of CO and horizontal gradients are 
significantly smaller, largely due to the coarser spatial resolution of MLS nadir measurements 
which were interpolated horizontally between orbital tracks. It is interesting to note that the 
large increase in CO in 2006 during October-November nearly disappears in December, yet O 3 
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in December in the western Pacific remains elevated (Figure 2c). Logan et al. [2008] attribute 
some of these increases in O3 due to lightning NOx in the Indonesian region. 

6 Impact of biomass burning during the 2006 El Nino 

In sections 4 and 5 we analyzed the El Nino related changes in tropospheric O3 and CO in 2006 
using both the model and observations. In this section we analyze the model results to estimate 
the relative importance of biomass burning and large scale transport in producing the El Nino 
related changes in tropospheric O3. Our approach is similar to the one used by Chandra et al. 
[2002] in analyzing the effect of the 1997 El Nino on tropospheric O 3 . The model was run in 
two modes: The first mode explicitly included the NOx and CO emission rates associated with 
the Indonesian fires as given in Table 1. The results of these runs are discussed in sections 4 and 
5. In the second mode, the model was run by excluding the contributions from the Indonesian 
fires. The model runs in this mode are used to assess the effects of large scale transport on 
tropospheric O3 as shown in the bottom panels of Figures 4a-4c. The difference of the two runs 
gives the contribution from biomass burning as shown in the upper panels of Figures 4a-4c. 

In October 2006 the increase in O3 due to biomass burning is about 4-6 DU (Figure 4a, upper 
panel). This increase is limited to a small region comprising the islands of Sumatra, Java and 
Borneo in the Indonesian region. It is comparable to 4-8 DU increase caused by changes in 
meteorological conditions (Figure 4b, lower panel). The latter, however, extends over a much 
larger area in the Indian Ocean encompassing regions of the Indian subcontinent and south-east 
Asia in the north, to western Australia in the south. These patterns are essentially similar in 
November 2006 (Figure 4b). In December, 2006, the O3 increase due to biomass burning is 
reduced to zero and the dynamical component of the O 3 field is identical to Figure 2c (lower 
panel). 

Table 2 shows the model simulations of monthly averaged global tropospheric burden of CO, O3 
and NOx with and without Indonesian fire emissions. The table shows that during October 2006 
the Indonesian fires contributed about 36 Tg which was about 10.7% of the global budget of CO. 
The increase in CO during November and December were respectively 12.6% and 8.1% of the 
global budget. These changes are significant and suggest that the El Nino related fires in the 
Indonesian region may have contributed significantly to global pollution. In comparison the 
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impact of fire was relatively small (< 3%) on the global O 3 burden. Logan et al. [2008] have 
suggested NOx production due to lightning as a contributing factor to the positive anomaly in 
TCO during December 2006. The NOx production due to lightning is implicitly included in the 
model calculations. In the tropical Pacific region (12°S-12°N, 90°E-170°E), lightning produced a 
total of 0.032 Tg NO in October 2006. This value increased to 0.0492 Tg NO in November and 
0.0562 Tg NO in December 2006. 

7 The El Nino of 2006: A long-term perspective 

Though the main focus of this study is the El Nino event of 2006, such events are very frequent 
and may have significant impact on global pollution. According to the World Meteorological 
Organization (WMO) an El Nino (La Nina) event is identified when sea surface temperature 
(SST) in the Nino 3.4 region (5°S-5°N, 120°W-170°W) is at least 0.5°C above (below) normal 
when averaged over three consecutive months. The relation between ASST and ATCO (the 
difference in TCO between eastern and western Pacific) is shown in Figure 5. The ATCO time 
series is derived by subtracting TCO in the eastern Pacific (15°S-15°N, 135°W-180°) from TCO 
in the western Pacific (15°S-15°N, 95°E-140°E). The time series is deseasonalized and smoothed 
with a three-month running average to be compatible with Nino 3.4 SST. TCO is derived from 
29 years (1979-2008) of total column ozone measurements from Nimbus 7 TOMS, Earth Probe 
TOMS and OMI using the convective cloud differential (CCD) method [e.g., Ziernke et al, 
2005]. Figure 5 shows a very robust correlation between ASST and ATCO. Even though the 
1997 El Nino produced largest perturbations in SST and TCO in the tropical Pacific, El Nino 
events like 2006 are more common. They are also episodic. Before 1997 they occurred about 
every 3-4 years, but more recently they have occurred every 2 years. A frequent occurrence of 
El Nino, similar to 2006, has the potential of increasing the global pollution triggered by the 
dryness and the resulting forest fires in the Indonesian region. 

8 Summary and Conclusion 

We have studied the effects of the 2006 El Nino on tropospheric O3, and CO at tropical and sub- 
tropical latitudes measured from the OMI and MLS instruments on the Aura satellite. The zonal 
characteristics of observed changes in these constituents are similar to those reported by Logan et 
al. [2008] based on measurements of these constituents from the TES instrument on the same 
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satellite. During October and December 2006, both these studies revealed a dipole like structure 
in TCO in the tropics with an increase over the western Pacific region and a decrease over the 
eastern Pacific region. The dipole structure was weaker during November 2006. 

The 2006 El Nino-induced drought allowed forest fires to spread rapidly during October and 
November in the Indonesian region. The effects of these fires are clearly seen in the 
enhancement of CO concentration measured from both the MLS and TES instruments. We have 
used a global model of atmospheric chemistry and transport (GMI CTM) to quantify the relative 
importance of biomass burning and large scale transport in producing observed changes in TCO. 
The model results show that during October and November both biomass burning and 
meteorological changes contributed almost equally to increases in TCO in the Indonesian region. 
The biomass component was 4-6 DU but it was limited to the Indonesian region where the fires 
were most intense. The dynamical component was 4-8 DU but it covered a much larger area in 
the western Pacific and Indian Ocean extending from South East Asia in the north to western 
Australia in the south. During December 2006 the effect of biomass burning was reduced to zero 
and the observed changes in TCO were mostly due to dynamical effects. Chandra et al. [2002] 
have obtained similar results for the 1997 El Nino using ozone measurements from Earth Probe 
TOMS and comparing them with the GEOS-CHEM model. The 1997 El Nino was significantly 
stronger than the 2006 El Nino and caused greater perturbations in TCO due to both forest fires 
and large scale transport. For the 2006 El Nino the model shows an increase of 2-3% in the 
global burden of tropospheric ozone. In comparison, the global burden of CO increased by 8 - 
12 %. 

This study shows that the 2006 El Nino can be characterized from long records of tropospheric 
O 3 and sea surface temperature as a moderate event, yet still the contribution to pollution from 
biomass burning emissions over Indonesia was substantial, both locally and also when evaluated 
on a global basis. In recent years El Nino events have occurred with greater frequency than in 
previous years dating back to 1979. Both the frequency and scale of emissions from biomass 
burning suggests that even relatively moderate El Nino events can be an important source of 
pollution in the troposphere. 
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latitude band 30°N-35°N show offset differences of 5-10 DU (with GMI larger than OMI/MLS). 
For the upper right panel in Figure A1 (corresponding to east China) the model and measured O 3 
in summer are comparable but the model values are larger than measured values by 5-10 DU in 
winter and spring months when STE is greatest. 


Martin et al. [2002] and Chandra et. at [2002] compared the seasonal variability in tropospheric 
O 3 from the GEOS-Chem model and TOMS CCD measurements. A large discrepancy was 
found between the model and the CCD measurements over low latitudes in northern Africa 
including Abidjan (5°N, 4°W). The model indicated maximum tropospheric column O 3 around 
December-F ebruary during peak time of the biomass burning whereas TOMS showed greatest 
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amounts around June-October at a time when there is little or no biomass burning in northern 
Africa. Our analysis shows that the OMI measurements agree much better with model than the 
TOMS measurements pertaining to this seasonal variability over Africa. The GMI model 
indicates that enhancements in tropospheric O 3 over northern Africa from the biomass burning 
lie almost entirely below 500 hPa and are highly localized to the relatively small region of 
burning emissions. 

Figure A2 shows statistical comparisons of model and observations. The three panels in this 
figure show their zonal correlations (top panel), mean zonal offset difference (middle panel) and 
mean zonal RMS (lower panel). The grey shading in the top panel represents areas that are not 
statistically significant at 99.9% confidence level. Figure A2 shows that the positive correlation 
between the model and OMI/MLS TCO, extending out to mid-latitudes, drops to negative values 
in higher latitudes. Several factors could contribute to near-zero or negative cross-correlations. 
Aside from known model uncertainties which include STE, OMI retrievals become less sensitive 
in measuring tropospheric O 3 in high latitudes (i.e., high solar zenith angles) and have reduced 
sensitivity in measuring O 3 in the lower troposphere, especially for high solar zenith angles and 
large slant columns. The model and observations compare the best in the tropics and subtropics 
extending out to latitudes of about ±30° to ±35°. This is the latitude range chosen for GMI and 
OMI/MLS measurements for studying the 2006 El Nino event. 
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Table 1: Emissions of CO and NO for 2006. 



CO (Tg) 

Indonesia 1 

World 1 

NO (Tg) 

Indonesia 11 

World’ 

July 

5.9 

92 

0.10 

7.3 

August 

16.6 

103 

0.29 

7.4 

September 

18.3 

99 

0.32 

6.8 

October 

44.2 

73 

0.76 

6.3 

November 

6.3 

79 

0.11 

6.3 

December 

0.2 

103 

0.00 

7.1 


d Biomass burning emissions only. 

^Emissions from all sources, including fossil fuels, biofuels, and biomass burning, except 
biomass burning from Indonesia. 
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'Fable 2. Monthly Averaged Global Tropospheric Burden of CO, O 3 , and NOx for the Model. 



Sep. 

Oct. 

Nov. 

Dec. 



CO (Tg) 


Indonesia 

351 

371 

371 

369 

No Indonesia 

336 

335 

329 

341 

ACO 

15 

36 

42 

28 

AGO (%) 

4.6% 

10.7% 

12.6% 

8.1% 



O 3 (Tg) 


Indonesia 

354 

351 

346 

332 

No Indonesia 

347 

343 

337 

324 

AO 3 

8 

8 

10 

8 

AO 3 (%) 

2 . 2 % 

2.4% 

2 . 8 % 

2.4% 



NOx (Tg N) 


Indonesia 

0.115 

0.121 

0.129 

0.133 

No Indonesia 

0.113 

0.117 

0.128 

0.133 

ANO x 

0.002 

0.004 

0.001 

- 0.000 

ANO x (%) 

2 .1% 

3.0% 

0 . 8 % 

- 0.3% 
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Simulations with (“Indonesia”) and without (“No Indonesia”) the 2006 Indonesia fire emissions. 


Figure Captions. 

Figure 1. (a) Tropospheric column ozone measured in DU from OMI/MLS for October 2004 
(top panel), October 2005 (middle panel), and October 2006 (bottom panel), (b) Same as (a) but 
instead for the GMI model, (c) The difference (GMI minus OMI/MLS) between Figure la and 
Figure lb. 

Figure 2. (a) Inter-annual difference (October 2006 minus October 2005) of tropospheric 

column ozone measured in DU from OMI/MLS (top panel) and the GMI model (bottom panel), 
(b) Same as (a) but instead for November, (c) Same as (a) but instead for December. 

Figure 3. (a) Inter-annual annual difference (October 2006 minus October 2005) of 215 hPa CO 
from MLS (top panel) and 227 hPa CO from the GMI model (bottom panel). All quantities 
represent volume mixing ratio in units ppbv. (b) Same as (a) but instead for November, (c) 
Same as (a) but instead for December. 

Figure 4. (a) Inter-annual difference (October 2006 minus October 2005) of GMI tropospheric 
column ozone (in DU) of the Indonesia biomass burning component (top panel) and the 
dynamically driven component (bottom panel), (b) Same as (a) but instead for November, (c) 
Same as (a) but instead for December. 

Figure 5. OMI/MLS time series ATCO (dark curve) plotted with Nino 3.4 time series ASST 
(light curve). ATCO is measured in DU while ASST is in Celsius degrees and was multiplied by 
a factor of three for scaling with ATCO. The ATCO time series is derived by subtracting TCO in 
the eastern Pacific (15°S-15°N, 135°W-180°) from TCO in the western Pacific (15°S-15°N, 95°E- 
140°E). The ATCO time series was deseasonalized and was also smoothed with a three-month 
running average to be compatible with the Nino 3.4 time series. 
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Figure Al. Monthly time series comparisons of GM1 (dotted curves) and OMI/MLS (solid 
curves) tropospheric column O 3 (in DU) for three separate regions: American (left four panels), 
African (middle four panels), and Indonesian (right four panels). 

Figure A2. (top) Zonal correlations between OMI/MLS and GMI tropospheric column O 3 . 
(middle) Calculated mean zonal offset difference (GMI minus OMI/MLS) measured in DU 
between model and observation, (bottom) Calculated zonal RMS of difference (GMI minus 
OMI/MLS) in DU. For both model and observation, fields for similar months (January - 
December) from 2004 to 2006 where averaged together to derive mean 12-month latitude- 
longitude gridded distributions of tropospheric O 3 . All statistical calculations were done zonally 
for each latitude band using the 72 gridded longitudinal values from model and observation. 
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